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MultiChann I Grating Design 

Field of the invention 

The present invention relates broadly to a multi-channel grating design method and to 
multi-channel grating structures fabricated utilising the multi-channel grating design method. 
Background of the invention 

Multi-channel grating structures are typically written into photosensitive waveguides. 
The grating structure comprises a refractive index profile induced in the photosensitive 
waveguide, which in turn determines the optical characteristics such as the reflection, 
transmission, and group delay characteristics of the resulting grating structure. 

The refractive index profile for a multi-channel grating structure is induced by applying 
a sampling function to a given single-channel grating profile, the single-channel grating profile 
typically being a periodic profile achieved by exposing the photosensitive waveguide to an 
appropriate light beam emerging from e.g. a suitable phasemask. 

hi a vast majority of previously reported work on multi-channel gratings, a so-called 
Sine-sampled design has been used. For the Sine-sampling approach, an TV-channel grating 
design can be obtained by a direct in-phase summation of N identical seeding gratings [with 
k(z) - amplitude grating amplitude, 6 (z) - grating phase] equally spaced in the frequency space: 

where 

Qsm «Xr =l oo ^-^-0Afe/23*^^sinc[^(Afc~2TO)/2],sinc(x) ^sm( x )fx, 
and Ak is the channel spacing. 

This design will be referred to as "in-phase" grating design herein after by the applicant. 

An example of this design is shown in Figures 9 (c) and (d) with corresponding spectral 
characteristics as shown in Figures 9 (a) and (b). The maximum value of the refractive index 
change required to implement this multi-channel grating design is given by a siitiple expression: 

An<r> = A^„ (2) 



08/24/2003 SUN 23:17 [TX/RX NO 5294] 



25. HUU. ZVUJ 10. \ J 



WO 02/069006 

PCT/AU02/00160 

2 

where Am, is the maximum refractive index change required for the single seeding grating. 
Since any photosensitive fiber used to fabricate Bragg gratings has material, limits of the 
maximum achievable photoinduced refractive index change An„ this represents a limitation on 
the maximum number of channels that can be recorded in a given fiber. Thus it is highly 
desirable to reduce a required An^ as much as possible. 

At least preferred embodiments of the present invention seek to provide an alternative 
multi-channel grating design in which the maximum refractive index change required as a 
function of the number of channels is reduced without deterioration of the grating spectral 
characteristics when compared with the prior art grating designs discussed above. 
Summary of the invention 

In accordance with a first aspect of the present invention there is provided a method of 
calculating a sampling function for fabricating a //-channel grating, the method comprising the 
steps of forming a summation of Aperiodic seeding functions each describing a refractive index 
variation, wherein each periodic function includes a phase shift value with respect to the other 
functions, and wherein at least one of the phase shift values is non-zero. 

Accordingly, the present invention can provide a "de-phased" grating design rather than 
the "in-phase" prior art grating design. 

Preferably, the summation of the N periodic functions comprises a Fourier analysis. The 
result of the Fourier analysis maybe expressed as: 

(-1 

where Q is the amplitude and y is the phase of the sampling function. 

The method may further comprise the step of deterrnioing a set of the phase shift values 
for which a maximum of the sampling function amplitude Q is minimised. 

Alternatively, the method further may comprise the step of deterrnining a set of the 
phase shift values for which a maximum difference between a minimum and a maximum of the 
sampling function amplitude Q is nnnimised. 
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Alternatively, the method may further comprise the step of determining a set of the 

i P hase shift values ft* which a mean-square-deviation in the sampling function : ampHturie Q is 

; minimised. 

j The step of determining the set of phase shift values may comprise direct scanning 

j through all combinations or conducting a variational analysis, or using other fonns of externum 

: search numerical techniques, or a simulated annealing Monte Carlo approach. 

The grating may be multi^dimensional, wherein the periodic seeding functions are multi- 
j dimensional. 

In accordance with a second aspect of the present invention, there is provided a method 
for fabricating a multi-channel grating comprising the step of calculating a sampling function in 
accordance with a method as defined in the first aspect. 

The multi-channel grating may e.g. be fabricated utilising photo-induced refractive index 
; changes in a photosensitive waveguide material, etching techniques, or epitaxial techniques, or 

| a developing technique such as a photo polymerisation process. 

I 

In accordance with a third aspect of the present invention, there is provided a multi- 
channel grating structure fabricated utilising a method of fabrication as defined in the second 

I aspect 

i 
i 
i 

Brief description of the drawings 

Preferred forms of the present invention will now be described with reference to the 

i accompanying drawings. 

j 

Figure 1 shows a four-channel grating design using the sampling function with 
minimised maximum amplitude embodying the present invention. 

Figure 2 shows the maximum value (maximum minimisation) of the sampling function 
amplitude for different number of channels, (a) Comparison with the worst case (in-phase 
sampling function design), (b) Comparison with the asymptotic theoretical value «Jn . 

; Figure 3 shows a 16-channel grating design using maximum minimisation approach 

j embodying the present invention. 

I 
i 
i 
i 

i 
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Figure 4 shows the maximum value (difference minimisation) of the sampling function 
amplitude for different number of channels. 

Figure 5 shows a 16-channel grating design using the difference minimization approach 
embodying the present invention. 

Figure 6 shows a 17-channel grating design using 1he variational minimization approach 
embodying the present invention. 

Figure 7 shows the maximum value (variational approach) of the sampling function 
amplitude for a different number of channels. 

Figure 8 shows an experimental set up for writing a multi-channel grating structure of a 
multi-channel grating design embodying the present invention. 

Figure 9 shows a prior art four-channel grating design using an in-phase-sampling 
function. 



Detailed description of the embodiments 

The preferred embodiment described provides a multi-channel grating design in which 
the maximum refractive index change is less than directly proportional to the number of 
channels N, and which avoids the presence of un-wanted side bands in the spectral 
characteristics of the resulting grating, thereby improving on prior ait multi-channel grating 
designs. 

In the preferred embodiment, a sampling function which periodically modulates the 
amplitude of a given single-channel grating (seeding grating) is utilised, similar to prior art 
multi-channel grating designs. However, in addition to the periodic modulation of the 
amplitude of the seeding grating, a periodic modulation of the phase of the seeding grating is 
also introduced. Accordingly, the resulting design function in the preferred embodiment may be 
expressed as: 

where the phase of the sampling function v^Hfa) and the sampling amplitude Q=Q(z) are 
given by: 



08/24/2003 SUN 23:17 [TX/RX NO 5294] 



PCT/AU02/00160 



e 2 (^) = 4^cos(a,^aJcos(n / Afe/2 + p i )cos(n i7 Afe/2 + P p ) > and 
(^(z) = tan H 



X)^ 2 sing, cosfo, Ate / 2 + fl) 
coscr, cos(n, Afe / 2 + ^ 



AT is even 



or 



(?(*) = <■ 2] cosa / C0S ( n / Afe/2+ ^)- 

1=1 



4 £cos(a, - ( ^^) c os(n 7 Z^fe/2 + p / )cos(» ;7 Afe/2+p ; ,) + l, 



and 



LSm"" 0 2oosa / cos(« f AAz/2 + p f +1) 
where n 2 ss2/-iV-l and s2^-J\f-l. 



Wis odd, 



In the above expressions for Q(z) and pffej we use notations 
a m d» r + )/2 7 /3 t m(fa- )/ 2 and set j* (Ar+1)ya = 0 for odd number of channels. Now for 
any given W there will be a set of {a,},{p f } (or equivalently a set of ) which optimizes the 
sampling function amplitude g(z) along the grating structure according to a selected 
optimization approach. By direct calculations it is straightforward to show that 



Q 2 dz~27iN/Ak, 



(4) 



for any choice of cti and ft. This expression, in turn, leads to an asymptotic formula for the 

minimum possible An N corresponding to an "ideal" situation when g(z)=V77 and only the 

grating phase is nontrivially modulated (by addition of an approprietly chosen phase of the 
sampling function iy(z) ): 

We note, that, in practice, the limit Q(z) - can be reached only approximately. For 
example, for the maximum minimization approach, one should solve a minimax problem and 
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find Q ma iz;a^\^)for which nuK,^<x;a™ P}^) }-«m to#lPl) a>K,fe(«:«/,Pf>}. 
To find the optimal set # for a relatively small number of channels one may use direct 

numerical scanning through all possible combinations of the dephasing angles. An example of 
the optimized "out-of-phase" design for a four-channel grating is shown in Fig. l(a>(d). For 
this design a, = 0.5759^ a 2 = Q 9 J3 X = J5 2 = 0 . This design reduces the maximum An N of the in- 
phase design by almost 40%. 

For N»l location of the minimizing set (aj opt> ,/}f 0pt) )is a nontrivial exercise. Even 
rough direct scanning through all possible sets of angles (followed by efficient numerical 
minimum search routines) quickly becomes numerically inefficient To solve the optimization 
problem for large N we use, in an example embodiment the so-called simulated annealing 
method - a Monte Carlo approach for minimization of multi-variable functions. This statistical 
method samples the search space in such a way that there is a high probability of finding an 
optimal or a near-optimal solution in a reasonable time. The term "simulated annealing" is 
derived from the analogy to physical process of heating and then slowly cooling a substance to 
obtain a crystalline structure. 

To start, the system state is initialized. A new configuration is constructed by imposing a 
random displacement If the energy of the new state is lower than that of the previous one, the 
change is accepted and the system is updated. If the energy is greater, the new configuration is 
accepted with some probability. This procedure allows the system to move consistently towards 
lower energy states, yet still jump out of local minimal due to the probabilistic acceptance of 
some upward moves. Results of the application of simulated annealing to the maximum 
minimization approach are summarised in Figure 2. 

Figure 2 shows that Eq. (5) gives a reasonable estimate for An N (20% accuracy or 

better) for N > 6. The growth factor to represents a significant improvement in 

comparison to linear ~7V dependence (almost a factor of 4 reduction for 16-channel gratings, 8 
for 64-chaimel ones and so on). An example of the optimal 16-channel grating design is shown 
in Fig. 3(a)-(d). 

Figures 1 and 3 demonstrate implementation of QJ^iz) reduction strategy without 
trying to avoid touching the zero level at some z. Zeros in the fibre Bragg grating (FBG) 
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amplitude may lead to the increased phase errors (appearance of phase jumps) and should be 



value of Q(z) along z, but in minimizing the maximum difference between the maximum and 
minimum of the sampling function amplitude. Mathematically this may be formulated as 



finding Q+tea^fijT 0 ) for which max,^(z;a^\pr ) )}-niin = fe rfto (z;a^°,p^>)}= 



This approach may be implemented by using the same simulated annealing algorithm, 
and the corresponding results are given in Fig. 4. An example of this version of optimal 
(difference minimisation) "out-of-phase" design for a 1 6-channel grating is shown in Fig. 5(a)- 
(d). The total minimization of A/x^ is slightly worse than in Fig. 3, but the zeros in the grating 
amplitude are avoided* 

It is noted that for both the- maximum minimisation and the difference minimisation 
approaches, points with the best optimisation quality may be used to get reasonable (but not 
exactly the best) optimisation for some higher channel numbers. For example, the point N=9 
with the particularly good optimisation quality provides us with an effortless optimisation for 
AM51 number of channels. 

Another embodiment will now be described, in which optimization by the functional 
minimisation (variational approach) is utilised. The key property of this embodiment is that it 
relies on estimate of some integral functional rather than time-consuming numerical scanning in 
z. Quantitatively, proximity of Q(z) to the theoretical limit -Jn can be characterised by mean- 
square-deviation, 



avoided. Thus, arguably, the best minimization strategy could be, not in reducing the maximum . 




(6) 



where (/(z))=2 / = Afc/27r£* f(z)dz Ideal optimisation of Q(z) corresponds to the 
achievement of the average Q « Vtf and the zero mean-square-deviatioji from this average 
value. Using expression (4) and assuming Q » Va^ , one can show that 



(7) 



where 
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For finding minima of AQ the most efficient strategy is again the use of the simulated 
annealing method. 

The major advantage of the optimisation based on the functional minimisation compared 
with direct scanning is the speed: integration over z is carried out analytically which saves lots 
of computer time. However, the quality of variational optimisation itself is usually not as good. 
That said, for odd number of channels variational optimisation leads to a sampling function 
without zeros in amplitude (similar to difference minimisation approach). An example for a 17- 
channel grating design (variational approach) is shown in Fig. 6(a)-(d). For even number of 
channels variational optimisation leads to a sampling function with zeros in amplitude (similar 
to maximum minimisation approach)- Fig. 7 shows, the maximum value (variational approach) 
of the sampling function amplitude for a different number of channels. 

The advantages of the dephasing optimisation approaches of the embodiments described 
include that they give a very significant reduction for the maximum required An# . Moreover, 
they are conceptually simple and relatively easy to obtain- Furthermore, the limitation of the 
number of dephased periodic functions to N for a JV-chaxmel grating design avoids the presence 
of un-desixed side bands, which can deteriorate the quality of spectral characteristics of the 
grating. The applicants have recognised that the limiting of the number of dephased periodic 
functions to N while allowing variations in the amplitude of the resulting sampling function can 
enable design of TV-channel gratings of improved quality when compared with prior art designs. 

The implementation of the multi-channel grating design of the preferred embodiment in 
a grating structure requires grating writing apparatus with high spatial resolution to be utilised. 
Therefore, in a grating writing apparatus relying on photoinduced refractive index changes, the 
apparatus preferably comprises a beam focusing means to reduce the size of the beam in the 
core of the photosensitive waveguide. 

Figure 7 shows an example experimental set up 50 for writing a multichannel grating 
52 into an optical fibre 54. The experimental set up 50 comprises an interferometer 56 which 
includes a first acousto-optic modulation 58 being operated Tinder an acousto-optic wave at a 
first frequency Q\ 9 as indicated by arrow 14, An incoming light beam 60 is incident on the first 
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acousto-optic modulator 58 at a first order Bragg angle. The operating conditions of the 
acousto-optic modulator 58 are chosen such that the. modulator 58 is under, driven, whereby, 
approximately 50% of the incoming beam 60 is diffracted into a first order beam 62, and 50% 
passing through the acousto-optic modulator 58 as an un-diffracted beam 64. The ua-dififcacted 
beam 64 is incident on a second acousto-optic modulator 66 of the interferometer 56 at a first 
order Bragg angle, whereas the beam 62 is not. Accordingly, the beam 62 passes through the 
second acousto-optic modulator 66 without any significant loss. 

The second acousto-optic modulator 66 is operated under an acousto-optic wave at a 
frequency Q2, which propagates in a direction opposing the direction of the acousto-optic wave 
in the first modulator 58 as indicated by arrow 68. After the second acousto-optic modulator 66 
the first order diffracted beam 70 and the beam 62 are frequency shifted in the same direction 
(e.g. higher frequency), but by different amounts i.e. €1\ v Q 2 » 

The beams 62, 70 are then brought to interfere utilising an optical lens 72, and the 
resulting interference pattern (at numeral 74) induces refractive index changes in the 
photosensitive optical fibre 54, whereby a refractive index profile, i.e. grating structure 52, is 
induced in the optical fibre 54, 

In Figure 7, the optical fibre 54 is translated along the interferometer at a speed v, as 
indicated by arrow 74. 

It will be appreciated by a person skilled in the art that the experimental set up 50 shown 
in Figure 5 can be utilised to write a multi-channel grating structure of a multi-channel grating 
design embodying the present invention through suitable control of the first and second acousto- 
optic modulators 58, 66, in conjunction with a suitable control of the speed v at which the 
optical fibre 54 is translated along the interferometer 56 at any particular time. The high spatial 
resolution required to implement the multi-channel design of the preferred embodiment is 
achieved in the set up shown in Figure 6 by utilising the optical lens 72, with the practical limit 
of the beam size in the focal plane preferably being of the order of the waveguide core size. 

It will be appreciated by a person skilled in the art that numerous variations and/or 
modifications may be made to the present invention as shown in the specific embodiments 
without departing from the spirit of scope of the .invention as broadly described. The present 
embodiments are, therefore, to be considered in all respects to be illustrative and not restrictive. 
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For example, multi-channel gratings can be fabricated on the basis of the multi-channel 
grating design of the present invention using various known grating /writing techniques, 
including one or more of the group of photo-induced refractive index variation in photo 
sensitive waveguide materials, etching techniques including etching techniques utilising a 
phasemask, and epitaxial techniques. Furthermore, while the preferred embodiments have been 
described in the context of 1 -dimensional Bragg gratings, the present invention does extend to 
multi-dimensional multi-channel gratings. Such gratings have applications e.g. as photonic 
bandgap structures. 
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